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Abstract 

The  research  presented  here  investigates  the  use  of  vibro-acoustic  methods  to  improve  the  performance  of  a  PEM  fuel  cell  by  enhancing  water 
removal  from  the  active  reaction  sites  within  the  fuel  cell.  Removing  the  water  increases  the  available  reaction  sites  and  thus  increases  the  available 
power  for  a  given  operating  condition.  To  examine  the  new  water  removal  methods,  first,  the  production  of  water  in  fuel  cells  and  current  water 
removal  methods  are  reviewed.  Then,  the  new  methods  are  proposed  that  are  based  on  structural  and  acoustical  excitation  of  the  stack.  Specifically, 
the  use  of  flexural  waves,  acoustic  waves  and  surface  waves  to  remove  water  from  a  fuel  cell  stack  are  examined.  Analytical  formulations  are 
given  in  order  to  calculate  the  excitation  frequency  and  amplitude  required  to  move  a  droplet  resting  on  a  vibrating  bipolar  plate.  Depending  on  the 
droplet  radius  and  other  parameters,  it  is  estimated  that  a  water  droplet  resting  on  a  bipolar  plate  can  be  moved  by  structural  displacement  levels 
as  low  as  1  pm.  The  different  approaches  to  droplet  removal  are  compared  in  terms  of  the  minimum  vibration  energy  required  per  droplet.  Water 
production  in  a  commercial  fuel  cell  stack  is  then  estimated  and  used  as  a  test  case  to  compare  the  power  required  to  effect  removal  of  a  certain 
number  of  droplets  with  the  amount  of  power  produced  by  the  stack.  It  is  shown  that  a  water  droplet  clogging  a  plate  channel  may  be  moved  with 
parasitic  power  requirements  as  low  as  21  mW.  For  each  method,  the  energy  required  to  effect  droplet  removal  is  quite  small,  although  among  the 
three,  the  use  of  surface  acoustic  waves  may  be  the  best  option  in  terms  of  minimal  vibration  energy  and  implementation  feasibility. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  capital  cost  associated  with  fuel  cell  (FC)  stacks  is  one 
of  the  major  roadblocks  to  FC  vehicle  commercialization.  The 
prices  of  prototype  FC  vehicles  currently  being  introduced  by 
various  automakers  are  in  the  vicinity  of  US$  1  million  [1],  a 
large  part  of  which  is  due  to  the  stack  component  and  manufac¬ 
turing  costs.  As  the  name  implies,  a  fuel  cell  stack  is  a  stack  of 
individual  fuel  cells  connected  in  series.  The  stacks  are  sized  to 
obtain  desired  current  and  voltage  (and  hence  power)  outputs. 
One  approach  to  reducing  stack  costs  is,  of  course,  to  reduce  the 
cost  of  the  stack  components,  some  of  which  are  made  with  pre¬ 
cious  metals.  Alternatively,  if  a  fuel  cell’s  performance  could  be 
improved,  the  stack  size  required  for  a  given  power  output  could 
be  reduced,  such  that  either  fewer  or  smaller  stack  components 
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could  be  used.  In  reality,  contributions  from  both  approaches 
will  most  likely  be  required.  The  focus  of  this  paper  is  on  the 
latter  approach:  improving  the  stack  performance  of  a  polymer 
electrolyte  membrane  (PEM)  fuel  cell  stack.  A  PEM  fuel  cell  is 
considered  in  this  work,  as  it  is  one  of  the  most  preferred  fuel 
cell  designs  in  industrial  use,  due  to  its  low  operating  tempera¬ 
ture  of  approximately  80  °C  [2].  The  performance  enhancement 
methods  presented  here  are  aimed  at  effectively  removing  the 
water  produced  as  a  result  of  chemical  reactions  in  a  PEM  fuel 
cell,  so  that  the  available  reaction  sites  will  be  increased.  In  the 
present  study,  the  water  of  interest  is  in  the  form  of  condensation 

[3]. 

Water  management  is  an  important  part  of  fuel  cell  opera¬ 
tion.  Water  is  the  product  of  the  electricity-producing  chemical 
reactions  inside  the  stack.  In  addition,  the  reactant  gasses  are 
often  pre-humidified  to  improve  their  permeability  through  the 
electrolyte  [4] .  As  a  result,  some  of  the  humidification  water  is 
carried  by  the  H+  ions  from  the  anode  to  the  cathode  through 
the  electrolyte  membrane  in  a  process  known  as  electro-osmotic 
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Fig.  1.  Sources  of  water  in  fuel  cells. 


drag.  Water  can  also  diffuse  from  the  cathode  to  the  anode  due 
to  concentration  and  pressure  gradients  across  the  membrane 
electrode  assembly  (MEA).  The  locations  and  motions  of  water 
in  a  single  fuel  cell  assembly  are  shown  in  Fig.  1. 

While  some  water  is  beneficial  to  the  stack  performance,  too 
much  water  can  lead  to  a  condition  known  as  “flooding.”  Flood¬ 
ing  in  fuel  cells  is  the  condition  when  water  clogs  the  channels 
of  the  bipolar  plates  and  reduces  the  flow  of  reactant  gases  avail¬ 
able  for  chemical  reactions  [5].  Moreover,  flooding  causes  some 
of  the  water  to  block  the  reaction  sites,  thus  affecting  the  rate  of 
chemical  reactions  and,  in  turn,  the  FC  efficiency. 

As  a  result  of  the  potential  problems  that  result  from  excess 
water,  various  researchers  have  proposed  methods  to  remove 
that  excess  water.  These  methods  are  based  on  water  removal  at 
either  the  anode  or  the  cathode.  A  review  of  some  of  these  meth¬ 
ods  is  given  in  References  [6,7] .  One  conventional  water  removal 
method  is  the  use  of  capillary  action  or  wicking,  in  the  gas  diffu¬ 
sion  layer  (GDF)  [8].  The  GDF  is  made  out  of  a  porous  material 
and  is  located  between  the  electrodes  (anode  and  cathode)  and 
the  bipolar  plates,  as  shown  in  Fig.  1 .  Water  is  drawn  through  the 
pores  in  this  material  by  capillary  action.  Another  method  for 
water  removal  is  to  supply  air  (or  oxygen)  to  the  cathode  at  a  rate 
significantly  greater  than  the  stoichiometric  ratio.  Typically,  the 
oxygen  excess  ratio  is  2.0;  that  is,  oxygen  is  supplied  at  a  rate 
twice  as  high  as  that  necessary  for  the  chemical  reactions.  One 
effect  of  the  higher- than- stoichiometric  ratio  is  that  the  excess 
gas  flow  can  be  used  to  “sweep  out”  the  excess  water  on  the 
cathode  side  of  the  cells  [9] .  Finally,  the  cell  electrodes  are  often 
coated  with  a  hydrophobic  material  in  order  to  reduce  the  elec¬ 
trodes’  affinity  for  water.  These  methods,  although  satisfactory, 
increase  the  fuel  cell  costs  in  terms  of  production  and  parasitic 
losses. 

Another  group  of  water  removal  methods  generally  referred 
to  as  “anode  water  removal”  approaches,  utilize  gradient  mech¬ 
anisms:  pressure  [10],  concentration  [11]  or  temperature  [12]. 
Water  is  driven  from  the  cathode  to  the  anode  side  of  a  given 
cell  by  creating  at  least  one  of  those  gradients  across  the  MEA 
and  bipolar  plates.  The  main  disadvantage  of  using  this  method 


is  its  detrimental  effect  on  the  ion  transportation  rate,  as  the 
water  is  diffusing  through  the  membrane  in  a  direction  opposite 
to  that  of  the  H+  ions  traveling  to  the  cathode.  Also,  the  excess 
water  accumulated  on  anode  side,  if  not  removed  effectively, 
may  reduce  the  amount  of  active  reaction  sites,  thus  affecting 
the  rate  of  reaction  and  again  reducing  the  system  performance. 

From  the  above  discussion,  it  is  evident  that  water  manage¬ 
ment  in  fuel  cells  is  an  important  engineering  topic.  The  field 
of  vibro-acoustics  may  offer  non-conventional  methods  to  real¬ 
ize  water  management  in  fuel  cells.  One  approach  to  realizing 
water  management  in  fuel  cells  that  has  already  been  investi¬ 
gated  is  the  use  of  vibro-acoustics  to  achieve  atomization  of 
water  droplets  in  the  stack  [13].  The  present  work  examines  an 
alternate  approach:  the  bulk  movement  of  water  droplets  within 
the  fuel  cells.  Through  vibro-acoustic  means,  it  may  be  possible 
to  move  water  droplets  from  the  MEA  to  other  locations  such 
as  drain  ports.  At  the  same  time,  this  may  be  accomplished  with 
no  moving  parts,  should  require  very  little  power  and  would 
require  minimal  or  even  no  changes  in  the  fuel  cell  designs 
to  enhance  water  movement  in  fuel  cells.  Furthermore,  these 
methods  may  prove  to  be  localizable,  in  that  they  can  be  imple¬ 
mented  as  needed  at  specific  locations  in  the  stack.  For  example, 
it  may  be  sufficient  to  “turn  on”  vibro-acoustic  actuators  on  only 
specific  locations  in  the  stack,  such  as  clogged  channels,  while 
actuators  at  non-problem  stack  locations  may  be  left  “off”  until 
needed.  In  that  fashion,  power  will  not  be  wasted  by  exciting 
non-problem  areas  of  the  stack.  In  contrast,  water  removal  meth¬ 
ods  such  as  the  use  of  significantly  higher  than  stoichiometric 
air  ratios  are  typically  applied  to  the  entire  stack,  rather  than 
specific  locations. 

The  next  section  presents  a  theoretical  background  for  three 
vibro-acoustic  methods  which  may  be  used  to  realize  water 
droplet  motion  in  a  FC  stack.  These  methods  may  potentially 
be  applied  to  both  the  cathode  and  anode  sides  of  individual  fuel 
cells.  The  formulations  are  aimed  at  estimating  the  amplitude 
of  vibration  and  consequently  the  minimum  vibration  energy 
required  to  move  a  droplet  from  one  location  to  another  in  a  bipo¬ 
lar  plate  channel.  The  three  methods  are  differentiated  in  terms 
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of  the  mode  of  excitation  used  and  include  (a)  generation  of  flex¬ 
ural  waves  in  the  bipolar  plate,  (b)  use  of  acoustic  waves  in  the 
flow  and  (c)  use  of  localized  surface  waves  in  the  bipolar  plates. 
A  general  description  for  implementing  each  of  these  methods 
is  given  in  the  respective  sub- section.  The  amount  of  energy  per 
droplet  is  then  calculated  for  each  approach.  An  approach  for 
estimating  the  number  of  water  droplets  that  must  be  removed  is 
then  presented.  That  approach,  although  somewhat  crude,  pro¬ 
vides  a  means  for  comparing  the  proposed  methods  in  terms  of 
the  minimum  overall  vibration  power  required  to  effect  water 
management  in  a  given  fuel  cell  stack  and  then  comparing  those 
values  to  the  overall  power  produced  by  the  fuel  cell.  A  1.2  kW 
fuel  cell  was  considered  for  this  comparison  study.  Finally,  to 
provide  perspective  on  the  potential  efficiency  improvements 
that  may  be  obtained  a  simple  method  to  estimate  the  effect  of 
flooding  on  overall  FC  performance  is  also  given. 

2.  Theoretical  background 

The  water  management  methods  proposed  in  this  work  are 
intended  to  effect  water  removal  by  generating  a  force,  Fdr0p, 
that  is  equal  to  the  force  of  adhesion  between  the  water  and  the 
supporting  surface,  F adh-  In  equation  form,  the  basis  for  water 
removal  may  be  expressed  as: 


contact  angle  cpo  and  the  radius  of  the  assumed  sphere,  R ,  both 
of  which  are  shown  in  Fig.  2.  The  contact  angle  is  a  property  of 
the  surface  on  which  the  droplet  rests.  The  bipolar  plate  is  gener¬ 
ally  coated  with  a  hydrophobic  material  such  as  Teflon  to  make 
water  removal  easier.  The  contact  angle  for  a  Teflonized  surface 
is  around  105°  [16].  The  percentage  of  the  spherical  volume 
that  makes  up  the  droplet  increases  with  the  droplet  angle.  For 
example,  the  extreme  case  of  cpo  =  1 80°  would  imply  a  perfect 
sphere  resting  on  top  of  a  surface.  The  effective  volume  Veff  and 
surface  area  Seff  of  the  spherical  droplet  cap  may  be  calculated 
using  the  relations  [17] 

f  nR?  \  9 

Veff  =  (  -y-  )  (i  -  COS  <por(2  +  COS  <po),  (3) 

and 

SQff  =  1txR2{\  —  cos^o)-  (4) 

In  this  work,  the  force  Fdr0p  required  to  induce  droplet  move¬ 
ment  may  be  generated  by: 

(i)  flexural  plate  waves; 

(ii)  acoustic  waves; 

(iii)  surface  waves. 


Fdrop  —  F’adh-  (1) 

The  force  of  adhesion  exerted  on  a  droplet  resting  on  a  surface 
is  given  by  [14,15] 

Fadh  =  CFwi  1  +  COS  (po)lV,  (2) 

where  aw  is  the  surface  tension  of  the  liquid  droplet  and  w  is  the 
width  of  the  droplet  base,  w  =  2R  sin(^o)-  Depending  on  the 
case  examined,  the  supporting  surface  may  be  a  membrane  or  a 
bipolar  plate.  Note  that,  in  order  to  initiate  droplet  movement, 
Fdrop  must  actually  be  greater  than  Fadh  •  However,  once  initiated, 
movement  at  a  constant  velocity  requires  only  that  Eq.  (1)  be 
satisfied  and  that  equality  will  be  taken  as  the  benchmark  for 
realizing  water  droplet  removal. 

A  droplet  resting  on  a  surface  is  shown  in  Fig.  2.  The  assump¬ 
tion  is  that  the  droplet  can  be  represented  as  the  intersection  of  a 
sphere  with  a  flat  surface.  The  amount  of  water  depends,  on  the 


Fig.  2.  Schematic  of  a  water  droplet  cap. 


The  three  methods  discussed  in  this  work  differ  in  the  way 
that  the  left  hand  side  of  Eq.  (1)  is  calculated. 

With  this  background  and  general  set  of  equations,  specific 
methods  for  water  removal  can  be  discussed.  It  should  be  noted 
here  that  some  of  the  methods  discussed  below  are  based  on  the 
works  of  other  researchers  who  have  investigated  the  phenomena 
of  droplet  movement  in  general  terms  or  for  other  applications. 
It  is  the  goal  of  the  present  research  to  apply  these  studies  to 
improve  the  product  water  removal  in  a  FC  system.  This  is 
done  by  tailoring  some  of  the  equations  derived  previously  and 
by  discussing  implementation  issues  with  regards  to  FC  sys¬ 
tems.  The  implementation  of  these  methods  may  include  issues 
such  as  the  type  of  actuator  used  and  its  location  in  a  FC  stack. 
Also,  the  methods  presented  may  potentially  be  used  on  both  the 
bipolar  plates  as  well  as  the  membrane  for  water  removal.  Of 
course,  additional  information  such  as  membrane  materials  and 
mechanical  properties  would  be  required  to  apply  these  meth¬ 
ods  to  different  components  of  the  fuel  cell  stack.  The  analytical 
models  presented  here  are  discussed  with  respect  to  parame¬ 
ters  such  as  droplet  volume  and  excitation  frequency  to  induce 
droplet  movement  on  a  bipolar  plate. 

2.7.  Droplet  movement  using  flexural  plate  waves  (FPW) 

The  phenomena  of  droplet  movement  on  plates  have  been 
presented  by  Scortesse  et  al.  [18]  and  Biwersi  et  al.  [19].  Those 
authors  examined  the  generation  of  flexural  waves  on  a  plate  by 
exciting  the  plate  at  frequencies  in  the  range  of  the  higher  flexural 
modes.  The  flexural  plate  waves  were  used  to  impart  energy  to 
resting  droplets  and  thus  cause  movement  of  the  droplets  in  the 
direction  of  the  traveling  wave,  as  illustrated  in  Fig.  3.  The  main 
forces  acting  on  a  droplet  due  to  flexural  wave  propagation  are 
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Fig.  3.  Droplet  movement  on  a  plate  by  exciting  flexural  waves. 


due  to  [20]  the  radiation  pressure 
Pv  =  -pco2u2  ^2  +  —  ^  sin 2(kx), 

the  capillary  pressure 

,  1  1 

Pc  —  &  \  — - h 


R 


X 


R 


y 


and  the  hydrostatic  pressure 


Fh  =  pgx. 


(5) 

(6) 

(7) 


where  k  =  2n!X  is  the  wave-number,  X  the  corresponding  wave 
length,  p  the  liquid  density  and  co  is  the  excitation  frequency  in 
rad  s_1 .  The  parameters  B  and  A  are  non-linear  effect  constants 
with  the  ratio  B/A  =  5.2  for  liquids,  u  the  vibration  amplitude 
in  meters  and  v  is  the  direction  of  propagation  of  the  flexural 
wave.  The  equation  of  capillary  pressure  is  derived  according  to 
Laplace  Law,  explained  in  Reference  [21].  In  this  equation,  a  is 
the  surface  density  of  the  liquid  droplet  and  Rx  and  Ry  are  the 
radii  of  curvature  on  the  horizontal  surface  of  the  plate  in  the 
respective  directions.  For  a  spherical  droplet,  Rx  =  Ry  =  R.  Note 
that  the  equation  for  radiation  pressure  is  time  independent  and 
second  order.  That  equation  is  obtained  by  relating  the  pres¬ 
sure  induced  by  the  flexural  waves  at  the  solid-liquid  interface, 
the  ambient  pressure  and  the  droplet  radius  of  curvature  [20]. 
The  harmonic  term  in  the  equation  indicates  the  assumption  of 
sinusoidal  motion  of  the  plate  in  flexure. 

Of  the  three  pressures  listed  above,  the  radiation  pressure 
is  observed  to  be  the  dominant  contributor  to  the  overall  forces 
exerted  on  a  resting  droplet  [22] .  Defining  the  pressure  amplitude 
P0r  required  to  move  a  water  droplet  as: 


1  ?2 

=  u 


The  force  Fx  exerted  on  the  droplet  by  radiation  pressure  can 
then  be  calculated  by  integrating  the  radiation  pressure  over  the 
angles  <p  and  6 ,  shown  in  Fig.  2,  as: 

rv  o 

Fx=  /  For  sin 2{kx)R2  sin 2  tp  cos#d#d<p.  (9) 

Jo  Jo 

Equating  this  force  to  the  adhesive  force  in  Eq.  (2)  and  solving 
for  the  pressure  amplitude  F0r, 


<7^(1  +  cos  cp)w 

fo°  fo 71  si n2(kx)R2  sin2  cp  cos  §  dd  dtp 


(10) 


Given  particular  droplet  parameters,  Eq.  (10)  can  be  solved  using 
the  Matlab®  function  ‘dblquad.’  The  resulting  F0r  can  then  be 
used  with  Eq.  (8)  to  solve  for  the  minimum  vibration  ampli¬ 
tude  required  to  generate  flexural  waves  that  will  cause  droplet 
movement  according  to: 


u  = 


(l/4)p(2  +  B/A)co2' 


(11) 


The  resultant  vibration  amplitudes  for  a  few  representative  cases 
are  presented  in  Section  3.1.  Once  the  vibration  amplitude  u  has 
been  obtained,  the  corresponding  vibration  energy  required  for 
droplet  movement  using  this  method  can  also  be  calculated.  This 
is  given  by  [23] 

Fyib  —  y^PujU  Veff,  (12) 

where  U  is  the  velocity  amplitude  of  vibration,  pw  the  water 
density  =  1000  kg  m-3  and  Veff  is  the  effective  volume  of  the 
droplet  given  by  Eq.  (3). 

A  number  of  means  are  available  for  exciting  a  bipolar 
plate  in  order  to  generate  the  flexural  waves  required  to  effect 
droplet  movement,  including  the  use  of  electrodynamic  shak¬ 
ers  or  piezoelectric  (piezo)  patches.  Each  method  has  its  own 
advantages  and  disadvantages:  a  shaker  is  easily  implemented, 
but  may  require  excess  space  that  is  not  available.  Piezo  patches 
are  quite  compact  methods  for  applying  forces  to  structures, 
but  must  be  appropriately  applied  and  attached  in  order  to  be 
effective.  One  other  technology  that  is  available  is  the  piezoce¬ 
ramic  inertial  actuator  [24].  That  device  is  essentially  a  miniature 
shaker  that  uses  piezoceramic  rings  to  generate  forces  between 
an  underlying  structure  and  an  inertial  reaction  mass. 

No  matter  which  method  of  excitation  is  used,  modal  analysis 
of  the  bipolar  plates  should  be  performed  to  determine  resonant 
frequencies  before  selecting  and  placing  actuators.  The  reason 
for  this  is  that  excitation  of  a  plate  at  a  resonant  frequency  is  gen¬ 
erally  much  more  effective  and  efficient  than  at  a  non-resonant 
frequency.  Exciting  a  lightly  damped  structure  at  resonance  will 
result  in  significant  displacements,  even  with  relatively  small 
force  levels.  In  addition,  modal  analysis  can  be  used  to  guide 
actuator  placement  and  to  avoid  the  problem  of  attaching  actu¬ 
ators  at  node  locations.  This  may  be  an  interesting  subject  to 
pursue  for  further  study. 

Regardless  of  which  excitation  source  is  used,  it  may  be 
attached  to  one  or  both  of  the  end  plates  in  a  FC  stack.  End  plates 
have  channels  cut  only  on  one  of  its  sides  and  they  enclose  all 
the  other  fuel  cells  in  the  stack.  It  should  be  ensured  that  the 
force  exerted  by  the  source  is  enough  to  achieve  the  required 
vibration  levels  for  droplet  movement  in  all  the  bipolar  plates  in 
the  stack. 


2.2.  Droplet  movement  using  modulated  ultrasound 

The  next  model  presented  is  based  on  the  analytical  formula¬ 
tion  given  by  Whitworth  et  al.  [25].  Those  authors  demonstrated 
the  movement  and  accumulation  of  polystyrene  particles  in 
water  in  the  presence  of  high  frequency,  high  intensity  acoustic 
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Fig.  5.  Droplet  movement  on  a  SAW  excited  plate. 


minimum  acoustic  pressure  required  to  move  a  droplet  may  be 
obtained.  This  acoustic  pressure  is  thus  given  by: 


Fig.  4.  Particle  movement  in  a  pseudo  standing  wave  field;  (a)  polystyrene  par¬ 
ticles  in  water  and  (b)  droplet  movement  in  air  in  a  bipolar  plate  channel. 


1  +  cos  cpo)w 


(15) 


waves.  The  authors  observed  that  that  movement  was  realized 
by  generating  a  pseudo  standing  (slowly  moving)  wave  field. 
A  schematic  of  the  experimental  setup  used  to  generate  such 
a  wave-field  is  provided  in  Fig.  4(a).  Two  piezo-actuators  are 
located  at  opposite  ends  of  a  channel  and  are  excited  at  slightly 
different  frequencies,  given  by/i  and/2,  where/i  >fi.  The  result¬ 
ing  acoustic  force  causes  movement  of  the  particles  in  the  host 
medium  from  the  low  frequency  side  to  the  high  frequency  side. 
The  maximum  speed  at  which  the  particles  move  is  given  by: 


CO  1  —  0)2 

k\  +  k2 


(13) 


where  cot  and  ki  refer  to  the  angular  frequency  and  wavenumber, 
respectively,  that  result  from  the  excitation.  The  acoustic  force 
acting  on  a  particle  or  droplet  in  the  presence  of  ultrasound  is 
given  by  [25] 


Fz  =  VP2 


Bv  + 


(14) 


where  P  is  the  acoustic  pressure  amplitude  and  V  is  the  droplet 
volume.  The  terms  ff*  and  /  are  the  adiabatic  compressibilities 
of  the  droplet  and  surrounding  medium,  respectively.  Here,  BY 
is  a  quantity  that  results  from  the  relative  motion  between  the 
polystyrene  particles  and  water.  The  value  of  Bv  is  estimated  by 
using  the  excitation  frequency  and  properties  of  the  particle  and 
the  host  such  as  densities  and  shear  viscosity  coefficients  [26] . 

To  implement  this  method  in  moving  a  droplet  resting  on  a 
bipolar  plate  channel,  consider  a  water  droplet  in  air,  instead 
of  polystyrene  particles  in  water.  Piezoelectric  speakers  may  be 
used  as  actuators  in  this  case  and  are  placed  at  the  ends  of  the 
channels  as  shown  in  Fig.  4(b).  The  force  Fz  given  by  Eq.  (14) 
and  shown  in  the  figure  is  required  to  overcome  the  adhesive 
forces  that  hold  the  droplet  at  a  particular  location.  It  is  important 
to  note  that  the  actuators  may  be  placed  such  that  they  do  not 
cover  the  entire  channel,  so  that  the  moving  droplets  will  have 
space  to  exit  the  channel.  Also,  this  will  ensure  that  the  actuators 
do  not  interfere  with  the  reactant  gas  flows. 

The  force  given  in  Eq.  (14)  is  responsible  for  the  movement 
of  a  particle  or  droplet  in  a  medium.  If  this  force  is  equated  to 
the  adhesive  forces  acting  on  a  droplet,  as  given  by  Eq.  (2),  the 


Once  the  acoustic  pressure  amplitude  is  known,  the  correspond¬ 
ing  velocity  amplitude  of  vibration  U  can  be  found  using  a 
plane-wave  assumption 

U=—,  (16) 

pc 

where  p  is  the  air  density  and  c  is  the  sound  speed  in  air.  The 
required  vibration  amplitude  is  estimated  by  taking  Eq.  (16)  and 
dividing  it  by  the  angular  frequency  corresponding  tof\.  These 
results  are  presented  in  Section  3.2.  Once  the  vibration  amplitude 
is  known,  the  required  vibration  energy  can  be  calculated  as  for 
the  case  of  the  standing  waves  using  Eq.  (12). 

2.3.  Droplet  movement  using  surface  acoustic  waves 
(SAWs) 

Surface  waves,  as  their  name  suggests,  are  high-frequency 
waves  localized  at  the  surface  of  the  material  in  which  they 
are  excited.  It  is  because  of  this  property  that  SAWs  may  be 
used  to  excite  and  move  a  small  water  droplet  resting  on  a  plate. 
There  are  various  types  of  surface  acoustic  waves  based  on  oper¬ 
ating  frequency,  depth  of  penetration  and  sensitivity  [27].  The 
most  widely  used  surface  wave-types  are  Rayleigh  waves,  Lamb 
waves,  Love  waves,  Shear  Horizontal  (SH)  waves  and  surface 
transfer  waves,  each  of  which  is  classified  according  to  the  spe¬ 
cific  type  of  acoustic  wave  used. 

Rayleigh  waves  are  one  of  the  most  widely  used  waves 
because  of  their  high  energy  concentration  on  the  surface  [28] 
and  better  entrainment  and  directivity  characteristics  [29] .  More¬ 
over,  Rayleigh  waves  are  considered  to  be  highly  sensitive  to 
surface  mass  loading  [30].  This  characteristic  may  be  advanta¬ 
geous  in  moving  small  volumes  of  liquid  resting  on  a  surface 
excited  by  Rayleigh  waves.  The  Rayleigh  SAW  radiates  a  wave 
component  into  a  liquid  droplet  when  the  SAW  propagation  sur¬ 
face  is  in  contact  with  that  liquid.  The  wave  component  is  normal 
to  the  plate  surface,  as  shown  in  Fig.  5.  That  figure  also  shows 
a  pair  of  inter-digital  transducers  (IDTs)  that  may  be  used  to 
excite  Rayleigh  waves.  IDTs  are  generally  used  to  excite  high- 
frequency,  high-intensity  surface  waves  [31].  As  a  first  step  to 
investigating  the  possibility  of  moving  a  droplet  using  Rayleigh 
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waves,  it  is  important  to  estimate  the  force  exerted  by  such  a 
wave  on  the  droplet  and/or  at  the  surface.  The  equation  of  radi¬ 
ation  pressure  generated  by  the  Rayleigh  waves  was  derived  by 
Gol’dberg  and  Naugol’nykh  [32].  The  basic  theory  and  wave 
equation  for  Rayleigh  waves  is  given  in  References  [33,34].  In 
one  of  his  studies,  Viktorov  derived  expressions  for  wave  veloc¬ 
ities  at  the  solid-liquid  interface  [35].  That  study  assumed  an 
infinite  layer  of  fluid  resting  on  the  excited  plate.  Furthermore, 
that  work  was  based  on  both  analytical  and  experimental  results. 
He  observed  that  in  case  of  Rayleigh  waves,  the  wave  velocities 
for  the  plate  in  the  transverse  direction  are  substantially  higher 
than  those  in  the  longitudinal  direction,  which  leads  to  higher 
energy  transfer  at  the  solid-liquid  interface  [36]. 

The  approach  used  in  the  present  study  will  involve  estimat¬ 
ing  the  vibration  amplitudes  required  to  move  a  droplet  resting 
on  a  SAW-excited  bipolar  plate.  To  do  this,  it  is  first  shown  that 
the  transverse  component  of  displacement  resulting  from  the 
excitation  of  a  plate  is  dominant  on  the  surface.  It  is  this  compo¬ 
nent  that  may  be  responsible  for  the  droplet  movement.  Then,  the 
equation  for  the  pressure  at  the  solid-liquid  interface  is  derived. 
Using  this  pressure  value,  the  corresponding  force  exerted  on 
a  liquid  drop  in  the  direction  normal  to  the  plate  surface  can  b 
estimated.  This  force  is  then  equated  to  the  adhesive  forces  given 
by  Eq.  (2)  in  order  to  estimate  the  vibration  amplitude  required 
to  move  a  droplet  resting  on  a  SAW  excited  plate. 

To  investigate  the  use  of  the  SAW  method  in  a  FC  system,  a 
bipolar  plate  is  considered  with  IDTs  mounted  at  the  ends.  The 
voltage  input  to  these  IDTs  will  result  in  a  proportional  force 
exerted  at  the  solid-liquid  interface.  Plane  harmonic  Rayleigh 
waves  are  shown  propagating  at  the  plate  surface  in  Fig.  5.  The 
components  of  longitudinal  and  transverse  displacements  in  the 
v  and  z  directions,  Ur  and  Wr  respectively,  are  given  by  [34,37] 


Ur  =  Ak R  e-^Rz  - 


R  e  5Rz  )  sin(kRV  —  cot), 


k2  4-  v2 
/cR  -t-  ar 


(17) 


and 


Wr  =  AqR  e-^Rz  - 


2£r 


k2  4-  v2 
kr  -t-  a  R 


e  SrZ  I  cos(Arv  —  cot),  (18) 


where  A  is  a  constant  proportional  to  vibration  amplitude;  Ar 
is  the  Rayleigh  wave  number  =  co/Vr;  co  is  the  excitation  angu¬ 
lar  frequency;  Vr  the  Rayleigh  wave  velocity;  qR  =  kR  —  A2; 
k\  =  co/V i;  sr  =  kR  —  A2;  At  =  co/Vt.  Here,  V\  and  Vt  are  the  phase 
velocities  of  longitudinal  and  transverse  waves  respectively.  V\ 
and  Vt  can  be  calculated  based  on  the  physical  properties  of  the 
material 


Vt  =  +7A  (19) 

and 

Vi  =  Vo  +  2  n)/p,  (20) 


where  p  is  the  material  density  and  /x  and  A  are  Lame  constants 
that  are  determined  by  the  material’s  Young’s  modulus,  E  and 


Poisson’s  ratio,  v,  according  to: 


E 

a  =  - , 

2(1  +  v)’ 


(21) 


and 


Ev 

(1  —  2v)(l  +  v) 


(22) 


In  the  present  study,  the  bipolar  plates  are  considered  to  be 
made  of  steel  with  properties  p  =  7700  kg  m-3  ,E  =  195  GPa  and 
v  =  0.28.  Using  those  values  with  Eqs.  (19)-(22)  results  in  the 
phase  velocities  in  longitudinal  and  transverse  direction  can  be 
calculated  to  be  V\  =  5.7e3  ms_1and  Vt  =  3.1e3ms~1,  respec¬ 
tively.  The  Rayleigh  wave  velocity,  Vr,  is  then  calculated  as  a 
function  of  transverse  wave  velocity  Vt  and  Poisson’s  ratio  v 
[38]: 


0.87  +  1.12v 


(23) 


The  magnitude  of  displacement  amplitudes  given  by  Eqs. 
(17)  and  (18)  are  plotted  as  a  function  of  the  normalized  depth 
in  Fig.  6  by  using  the  parameters  given  above.  This  plot  shows 
that  the  Rayleigh  waves  are  concentrated  at  the  surface  and  that 
they  decay  in  the  material  as  a  function  of  the  wavelength.  The 
magnitudes  are  normalized  with  respect  to  the  amplitudes  at 
the  surface  (z  =  0)  and  the  depth  is  normalized  with  respect  to 
Rayleigh  wavelength  Ar.  Note  that  the  transverse  component 
of  displacement,  Wr,  is  large  at  the  surface  as  compared  to  the 
longitudinal  component,  Ur.  As  mentioned  earlier,  this  char¬ 
acteristic  may  prove  to  be  advantageous  in  exerting  high  force 
magnitudes  at  the  interface.  To  estimate  the  force  exerted  on  the 
droplet  resting  on  the  bipolar  plate  surface,  the  pressure  gener¬ 
ated  by  the  Rayleigh  waves  is  calculated  first. 

To  calculate  the  pressure  at  the  interface  (z  =  0),  the  pressure 
release  boundary  condition  may  be  applied,  which  is  given  by: 


(24) 
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where  ozz  is  the  normal  stress  component  in  the  z  direction,  given 
by: 


azz  —  ^ 


320  320 

dx 2  +  3  z2 


0  /  320  3  V 

+  dz2  dxdz 


(25) 


where  \fr  and  <p  are  the  scalar  and  vector  potentials,  respectively 
at  the  interface  given  by: 

(p  =  —A  ?  (26) 

and 

f  =  iA  ^RgR  ei(kRx-cot)-sRz  _  (27) 

4+4 

Substituting  Eqs.  (25),  (26)  and  (27)  in  Eq.  (24),  the  magnitude 
of  the  pressure  at  the  interface,  Pmt,  is  obtained  as 

D  4  — 2Ak2kR+AA2k2— 4/xkR  +  2/xkRk2+4yuARgRSR  2 

Pint  =  A  ot2  _  j.2  77'r  ' 

Z,/CR  /ct 

(28) 


Fig.  7.  Droplet  movement  requirements  on  vibration  amplitude  using  flexural 
plate  waves. 


The  software  package  Maple®  was  used  to  evaluate  Eqs. 
(24)-(28).  The  corresponding  force  exerted  by  this  pressure  at 
the  interface  can  then  be  calculated  by  multiplying  Eq.  (28)  by 
the  contact  area  of  the  droplet.  This  force  which  is  responsible 
for  droplet  movement  was  then  equated  to  the  adhesive  forces  on 
the  droplet  in  Eq.  (2)  and  solved  for  the  constant  A.  The  value  of  A 
was  subsequently  substituted  into  Eq.  (18)  to  estimate  the  trans¬ 
verse  vibration  amplitudes  required  to  move  a  droplet  resting 
on  a  SAW-excited  plate.  The  results  for  the  estimated  vibration 
amplitude  are  presented  in  a  later  section  of  this  work.  Once  the 
transverse  vibration  amplitudes  are  known,  the  corresponding 
vibration  energy  required  for  droplet  movement  can  again  be 
calculated  using  Eq.  (12). 

Note  that  the  techniques  discussed  above  are  for  moving  a 
single  droplet  of  water.  In  an  actual  fuel  cell  stack,  there  may 
be  several  such  droplets.  A  later  section  of  this  work  discusses 
the  method  and  assumptions  made  for  estimating  the  number  of 
droplets  in  a  particular  fuel  cell.  Once  the  number  of  droplets 
produced  per  minute  is  estimated,  the  three  proposed  methods 
are  compared  in  terms  of  water  removed  and  minimum  vibration 
energy  required  to  remove  this  water.  Only  water  droplets  at  the 
cathode  are  considered  in  this  study. 

3.  Results  and  discussion 

The  results  for  estimated  vibration  levels  for  each  of  the 
methods  discussed  above  are  presented  in  the  following  three 
sub- sections.  The  results  are  presented  in  terms  of  vibration 
amplitudes  required  to  move  a  single  droplet  clogging  or  block¬ 
ing  the  channels  of  a  bipolar  plate.  An  underlying  assumption 
here  is  that  the  water  droplets  form  in  the  channels  and  there  is  no 
water  flow  through  the  gas  diffusion  layer.  This  means  that,  for 
the  purposes  of  this  work,  FC  performance  is  characterized  by 
the  water  flow  in  the  channels.  Also,  it  should  be  noted  that  the 
water  droplets  will  continue  to  move  across  the  plate  surface, 
so  long  as  the  bipolar  plate  is  excited.  An  estimate  of  droplet 
production  in  a  specific  fuel  cell  system  is  presented  in  Section 


3.4  and  is  followed  by  the  calculation  of  the  minimum  vibra¬ 
tion  power  required  to  effect  droplet  movement  for  that  level  of 
droplet  production. 

3.1.  Method  1:  water  removal  by  flexural  plate  wave 
excitation  (FPW) 

The  estimated  flexural  plate  wave  vibration  amplitudes 
required  to  move  a  droplet  on  a  vibrating  plate  are  shown 
as  a  function  of  droplet  size  in  Fig.  7  for  three  differ¬ 
ent  excitation  frequencies.  The  following  numerical  values 
are  considered  in  plotting  the  results:  speed  of  sound  in 
water  =  1481  m  s-1 ;  vo  (see  Fig.  2)  =  2.5  mm;  surface  tension  for 
water-air  system  [8]  =  72  mJ  m-2;  B/A  =  5.2  for  liquids;  density 
of  water  =  p  =  1000  kg  m-3 . 

Note  first,  from  Fig.  7  that  while  the  vibration  amplitudes 
required  to  move  the  droplet  are  not  substantial,  the  correspond¬ 
ing  acceleration  levels  required  will  be  sizeable  because  of  the 
high  excitation  frequencies.  An  approximation  of  the  accelera¬ 
tion  levels  may  be  obtained  by  squaring  the  excitation  frequency 
(in  rads-1)  and  multiplying  it  by  the  displacement  amplitude. 
For  example,  given  a  2  mm  droplet  radius,  a  vibration  amplitude 
of  150  [xm  is  required  to  effect  droplet  movement  for  an  exci¬ 
tation  frequency  of  10  kHz.  That  corresponds  to  an  acceleration 
amplitude  of  almost  60,000  g.  Whether  these  high  accelera¬ 
tion  levels  can  be  attained  in  a  FC  bipolar  plate  is  a  subject 
of  future  study.  Specifically,  the  stresses  in  the  plates  must 
be  estimated  and  compared  to  the  yield  stress  for  the  partic¬ 
ular  plate  material  under  consideration.  Issues  associated  with 
fatigue  must  also  be  considered.  Furthermore,  it  is  possible  that 
such  high  acceleration  levels  or  even  displacement  levels  may 
damage  and/or  delaminate  the  catalyst  layer.  Such  problems  will 
require  consideration  in  future  works.  Even  more  important, 
experimental  verification  of  these  predictions  is  needed  before 
it  can  be  determined  without  a  doubt  that  such  an  approach 
is  even  feasible.  For  the  moment,  however,  the  real  point  of 
interest  is  in  determining  the  vibrational  energy  required  to 


V.  Palan  et  al.  /  Journal  of  Power  Sources  161  (2006)  1116-1125 


1123 


Fig.  8.  Effect  of  droplet  radius  on  vibration  amplitude  with  the  use  of  modulated 
ultrasound.  Frequency  difference  between  the  two  actuators  is  200  Hz. 

effect  droplet  removal  using  this  technique  as  will  be  presented 
shortly. 

3.2.  Method  2:  droplet  movement  using  modulated 
ultrasound 

Fig.  8  shows  the  results  obtained  for  required  vibration  ampli¬ 
tude  to  move  a  droplet  using  modulated  ultrasound.  The  refer¬ 
ence  actuator  was  excited  at  a  frequency  (f\)  of  3  MHz.  For 
this  method  to  work,  the  excitation  frequencies  have  to  be  very 
high.  At  lower  frequencies,  it  was  observed  that  substantially 
higher  vibration  amplitudes  are  required.  In  plotting  various 
results  in  the  figure,  the  following  numerical  values  were  con¬ 
sidered:  frequency  difference  f\  —  /2  =  200Hz;  reference  pres¬ 
sure  JPref  =  20e— 6  Pa;  contact  angle  cpo  =  105°;  surface  tension 
for  water-air  system  a  =  72mJm-1;  adiabatic  compressibility 
for  water  /3*  =  4.6e— 10  Pa-1;  adiabatic  compressibility  for  air 
/3  =  7.02e— 6Pa  ;  density  of  water  p  =  1000  kg  m  ,  density 
of  air  p=  1 .21  kg  m-3 ;  sound  speed  in  air  c  =  343  ms-1.  The 
reason  for  selecting  a  frequency  difference  of  200  Hz  was  that 
the  effect  of  frequency  difference  on  vibration  amplitude  was 
found  to  be  insignificant  above  200  Hz. 

As  shown  in  Fig.  8,  the  required  vibration  amplitudes  are 
very  low  as  compared  to  the  first  method.  It  is  clear  from  these 
observations  that  the  applicability  of  this  method  would  greatly 
depend  on  the  ability  and  feasibility  of  inducing  such  high- 
frequency  sound  into  the  fuel  cell  system.  Further  work  in  the 
actuator  and  control  mechanisms  required  to  induce  high  fre¬ 
quency  acoustic  waves  in  plate  channels  may  give  additional 
insight  into  the  viability  of  this  method.  Also,  as  with  the  first 
method,  the  analysis  conducted  here  is  also  for  moving  a  single 
drop. 

3.3.  Method  3:  droplet  movement  using  surface  acoustic 
waves  (SAW) 

The  vibration  levels  required  for  the  use  of  surface  acoustic 
waves  are  shown  in  Fig.  9.  An  excitation  frequency  of  10  kHz  is 


Fig.  9.  Droplet  movement  requirements  for  vibration  amplitude  using  Rayleigh 
waves  approach. 

considered  for  this  particular  case.  Note  that  the  vibration  levels 
required  for  droplet  movement  using  this  method  are  much  lower 
than  those  required  for  the  use  of  flexural  waves,  while  they  are 
larger  than  those  required  for  the  use  of  ultra-sound  waves.  One 
of  the  reasons  for  smaller  required  displacements  of  the  SAW- 
based  approach,  as  compared  to  the  flexural  wave  approach,  is 
that  the  surface  acoustic  waves  are  localized  to  the  surface  of  the 
solid,  while  in  the  case  of  the  flexural  waves,  the  entire  structure 
is  in  motion. 

3.4.  Calculation  of  water  droplet  production 

As  mentioned  in  Section  2.1,  the  three  methods  for  droplet 
movement  presented  here  can  be  compared  based  on  the  mini¬ 
mum  vibration  power  required  to  move  a  certain  number  of  water 
drops  per  unit  time.  In  order  to  calculate  the  vibration  power 
levels,  therefore,  some  information  on  the  amount  of  water  pro¬ 
duction  in  a  fuel  cell  is  necessary.  To  that  end,  water  production 
in  a  commercially  available  1 .2  kW  fuel  cell  system  was  consid¬ 
ered.  Although  a  chemical  reaction  balance  and  analysis  would 
give  the  exact  amount  of  water  produced  for  this  fuel  cell  [39], 
such  a  balance  was  not  done  here.  Instead,  data  from  the  user’s 
manual  of  the  fuel  cell  system  was  used.  According  to  that  doc¬ 
ument,  the  water  production  on  the  cathode  side  of  the  bipolar 
plates  in  the  fuel  cell  is  14.5  ml  min-1.  For  the  purposes  of  the 
present  discussion,  it  is  assumed  that  water  is  produced  and 
moves  in  droplet  form,  as  opposed  to  continuous  flow  of  water. 
Assuming  2  mm  droplet  radii,  such  that  the  effective  volume  of 
each  droplet  is  0.023 1  ml,  the  number  of  droplets  forming  on  the 
surface  of  the  bipolar  plates  per  minute,  Adr0p,  is  628.  That  is, 
during  every  minute  of  operation,  there  are  approximately  628 
drops  of  water  that  may  be  responsible  for  flooding  or  block¬ 
age  in  the  fuel  cell  stack  under  consideration.  It  is  reasonable 
to  expect  that  most  of  that  water  would  naturally  drain  from 
the  fuel  cell.  However,  to  provide  a  conservative  estimate  of  the 
power  required  to  effect  droplet  removal,  it  is  assumed  that  all 
of  that  water  will  remain  inside  of  the  stack.  Using  that  number 
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and  the  minimum  vibration  energy  required  to  move  a  droplet 
by  the  three  methods  discussed  earlier,  the  total  vibration  energy 
required  to  move  all  the  droplets  can  be  calculated.  This  com¬ 
parison  and  analysis  is  presented  in  the  examples  discussed  in 
the  next  section. 

As  mentioned  earlier,  the  assumptions  made  in  the  estimation 
of  number  of  droplets  provide  a  sound  platform  for  comparison 
of  the  above  discussed  water  removal  methods.  It  is  important 
to  remember  that,  as  discussed  in  Section  1,  the  goal  of  realizing 
effective  water  removal  is  to  improve  the  overall  FC  perfor¬ 
mance.  The  easiest  method  of  quantifying  potential  performance 
improvements  is  to  compare  the  predicted  stack  current  out¬ 
put  increases  achieved  when  implementing  each  of  the  different 
water  removal  methods.  The  increased  stack  power  can  then  be 
compared  to  the  power  required  to  effect  the  water  removal.  In 
the  end,  an  experimental  study  will  be  required  to  truly  vali¬ 
date  the  results.  Such  an  approach  is  under  consideration  for 
a  future  work.  At  this  preliminary  stage,  however,  in  order  to 
determine  if  such  an  experimental  study  is  warranted,  a  sim¬ 
pler  approach  will  be  taken.  The  assumption  is  made  that  the 
performance  of  a  FC  system  is  proportional  to  the  total  surface 
area  of  the  active  reaction  sites  [40] .  It  should  be  noted  that  this 
assumption  was  made  to  simplify  the  discussion  and  the  result¬ 
ing  calculations.  A  more  realistic  assumption  may  be  to  consider 
just  the  channel  area  as  opposed  to  the  entire  active  reaction 
area.  Such  assumptions  will  be  considered  in  future  works.  The 
fuel  cell  stack  under  consideration  has  50  fuel  cells  connected 
in  series  and  hence  50  active  reaction  areas.  Considering  each 
active  area  to  be  50  mm  x  50  mm,  the  total  active  reaction  area 
will  be  125e— 3  m2.  The  628  droplets  of  2mm  radius  respon¬ 
sible  for  flooding  would  thus  reduce  the  active  reaction  area 
by  7.9e— 3  m2.  Therefore,  the  percent  reduction,  Pred,  in  per¬ 
formance  may  be  estimated  as  Fred  =  (7.9/125)  x  100  =  6.3%. 
Note  that  this  is  just  the  effect  on  performance  every  minute.  If 
excess  water  accumulates  over  time,  Pred  may  increase  dramat¬ 
ically.  Although  it  is  understood  that  this  is  a  rough  estimation, 
the  6.3%  reduction  in  performance  demonstrates  the  effect  that 
excess  water  can  have  on  overall  FC  performance.  Furthermore, 
that  value  provides  a  benchmark  for  comparing  the  different 
water  removal  methods;  in  order  for  any  one  of  the  methods 
to  be  a  viable  candidate  for  implementation  in  a  fuel  cell,  the 
respective  decrease  in  Pre d,  must  outweigh  the  power  require¬ 
ments  of  that  particular  method. 

3.5.  Calculation  of  power  requirements 

The  assumption  is  that  628  drops  with  2  mm  radii  are  formed 
every  minute  or  a  2  mm  droplet  production  of  approximately 
10.5  drops  s-1.  Fig.  10  shows  the  minimum  vibration  energies 
required  for  the  three  proposed  methods  for  different  droplet 
radii.  The  variation  in  vibration  levels  for  the  three  methods 
arises  from  the  inherent  difference  in  the  way  motion  is  generated 
by  these  methods.  Methods  1-3,  as  discussed  earlier,  use  flexu¬ 
ral  waves,  acoustic  waves  and  surface  waves,  respectively.  Note 
that  even  though  the  required  vibration  amplitudes  for  droplet 
movement  were  lowest  while  using  modulated  ultrasound,  the 
vibration  energies  are  lowest  for  the  SAW-based  method.  This  is 


10 


10  r 


^  10% 
CD 
<D 

5  Iff1 

c 

O 

TO 

~  io*2 


10 


10 


-3 


Flexural  waves  @10kHz 
Acoustic  waves  @3MHz 
Surface  waves  @10kHz 


1.5  2  2.5  3  3.5  4  4.5  5  5.5 

Droplet  Radius  (mm) 

Fig.  10.  Minimum  vibration  energy  required  to  move  628  water  droplets  on  a 
bipolar  plate  using  flexural  waves,  acoustic  wave  and  surface  waves. 


because  of  the  high  frequencies  considered  in  the  former  case. 
From  the  figure,  it  appears  that,  theoretically,  the  use  of  surface 
waves  may  be  the  best  option.  Whether  it  is  feasible  to  gen¬ 
erate  such  waves  in  an  actual  fuel  cell  is  a  subject  of  further 
investigation. 

To  examine  the  overall  effect  of  the  droplet  removal  meth¬ 
ods  on  the  fuel  cell  operation,  the  power  to  implement  each 
method  can  be  determined  using  the  approximate  droplet  rate 
of  10.5  drops  s-1.  For  example,  using  flexural  waves  requires 
approximately  22  mJ  drop-1  or  roughly  230  mW.  Assuming  a 

2  mm  droplet  radius,  the  power  requirements  for  the  acoustic 
and  surface  wave  approaches  are  21  and  0.063  mW,  respec¬ 
tively.  Those  power  requirements  should  be  compared  against 
the  rough  estimate  of  a  drop  of  approximately  6.3%  or  75  W  in 
a  1200  W  system.  Clearly,  on  paper,  each  of  the  three  methods 
indicates  a  great  deal  of  potential  for  improving  the  overall  fuel 
cell  system  performance,  with  the  surface  acoustic  waves  offer¬ 
ing  the  greatest  potential  gain.  However,  for  each  method,  the 
power  requirements  are  so  small  in  comparison  to  the  poten¬ 
tial  power  improvements  so  as  to  be  negligible.  In  that  case,  the 
choice  of  which  water  removal  method  to  use  should  be  based 
on  implementation  feasibility  and  the  cost  of  the  added  sen¬ 
sors  and  actuators.  As  noted  before,  at  this  point,  the  analysis  is 
entirely  theoretical.  The  next  step  in  this  area  is  to  implement 
different  water  removal  methods  on  a  functioning  fuel  cell  and 
to  obtain  experimental  validation  of  the  predicted  performance 
improvements. 

4.  Summary 

The  present  study  was  motivated  by  the  fact  that  fuel  cell 
technology  is  still  in  its  embryonic  stages  when  it  comes  to 
its  broad-based  commercialization.  The  issue  of  performance 
enhancement  is  a  subject  of  perennial  research;  hence  the  focal 
point  of  the  present  study.  Higher  performance  results  in  smaller 
stack  size  and  hence  lower  costs.  Effective  water  removal  is 
not  only  essential  to  the  smooth  running  of  fuel  cells,  but  also 
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improves  the  performance  if  done  effectively.  To  that  end,  funda¬ 
mental  methods  were  investigated  that  use  waves  such  as  flexural 
plate  waves  (FPW),  acoustic  waves  and  surface  waves  (SAW) 
to  enhance  the  water  removal  process.  Analytical  models  were 
developed  to  estimate  the  vibration  levels  and  minimum  vibra¬ 
tion  energies  required  at  a  particular  frequency  to  move  a  water 
droplet  of  a  particular  size  and  volume  on  the  surface  of  the  bipo¬ 
lar  plate.  Specifically,  it  was  shown  that  droplet  movement  can  be 
induced  by  displacement  amplitudes  that  are  less  than  or  equal 
to  1  p,m.  Moreover,  it  was  found  that  excitation  by  surface  waves 
provides  the  best  option  in  terms  of  vibration  energies  required 
to  move  a  droplet  on  the  bipolar  plate.  In  summary,  it  was  the¬ 
oretically  shown  that  the  proposed  approaches,  if  implemented, 
may  potentially  result  in  better  performance,  hence  reducing  the 
overall  FC  stack  cost.  In  the  future,  experimental  or  finite  ele¬ 
ment  modeling  may  be  conducted  to  verify  whether  the  required 
displacement  and  acceleration  levels  can  be  obtained  in  a  fuel 
cell  stack.  Since  the  membranes  between  the  plates  are  con¬ 
structed  out  of  a  relatively  soft  material,  it  seems  plausible  that 
such  vibration  levels  can  be  achieved. 
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